Huge excitonic effects in layered hexagonal boron nitride 
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The calculated quasiparticle band structure of bulk hexagonal boron nitride using the all-electron 
GW approximation shows that this compound is an indirect-band-gap semiconductor. The solution 
of the Bethe-Salpeter equation for the electron-hole two-particle Green function has been used to 
compute its optical spectra and the results are found in excellent agreement with available exper- 
imental data. A detailed analysis is made for the excitonic structures within the band gap and 
found that the excitons belong to the Frenkel class and are tightly confined within the layers. The 
calculated exciton binding energy is much larger than that obtained by Watanabe et aZQj using a 
Wannier model to interpret their experimental results and assuming that h-BN is a direct-band-gap 
semiconductor. 

PACS numbers: 71.15.Mb, 71.35.Cc, 78.20.-e 
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Hexagonal boron nitride (h-BN) is one of the most 
anisotropic layer compounds and represents an inter- 
esting quasi-two-dimensional insulator analog to semi- 
metallic graphite. Due to its high thermal stability, 
h-BN is a widely used material in vacuum technol- 
ogy. It has been employed for microelectronic devices 0, 
for x-ray lithography masks 0, and as a wear-resistant 
lubricant 1^1 . The interest in h-BN has been renewed by 
the possibility of preparing boron nitride nanotubes that 
are far more resistant to oxidation than carbon nanotubes 
and therefore suited to high temperature applications. 
As their band gaps are predicted to be weakly dependent 
on helicity and tube diameter 0.0]. unlike tubes made of 
carbon, it has the potential of revolutionizing the elec- 
tronics industry. To understand the quasiparticle prop- 
erties of BN nanotubes one has first to have a complete 
understanding of the electronic and optical properties of 
bulk h-BN since BN nanotubes can be viewed as a rolled 
up BN sheets. 

Despite the large number of experiments pi IM Ifl ITo| 
devoted to the study of the electronic properties of bulk 
h-BN, both the direct and the indirect band-gap are not 
yet accurately known, i.e., values ranging from 3.2 to 
5.97 eV have been reported in the literature. In par- 
ticular, the latest results were obtained by Watanabe et 
who showed that this material is very promising 
for ultraviolet laser devices and assumed, based on the 
strong luminescence peak observed at 5.765 eV, that it is 
a direct-band-gap semiconductor. In addition, they in- 
ferred a gap of 5.971 eV in contradiction with our results 
and other GW quasiparticle calculations [Til Il2l Il3 | . 

In this Letter, we use our all-electron GW 
appr oximat ion 1 1 4j to study for the first time the opti- 
cal properties of bulk h-BN including electron-hole inter- 
action effects and show that the previous experimental 
interpretation of the excitonic structure in the band gap 
is not correct. To perform this study we used state of the 
art methods, where both the self-energy effects and the 



solution of the Bethe-Salpeter equation^] to compute 
the optical spectra from first principles [lfj. In particu- 
lar, our study confirms that this material is an indirect- 
band-gap semiconductor in agreement with previous GW 
calculation [TIL [T3. IT^ | . and that the interpretation of the 
optical spectra of h-BN, upon which the band gap of 
Watanabe et aZp] is based, is very intricate as they are 
dominated by strong electron-hole interaction effects. We 
show that theoretical optical spectra are reliable and are 
in excellent agreement with experiment |10|. In addition, 
we show that the exciton in the band gap belongs to a 
Frenkel class instead of a Wannier class as proposed by 
Watanabe et al 1,] from the analysis of their experimen- 
tal data, and we found that its binding energy is huge 
compared to theirs. 

The crystalline structure of h-BN is hexagonal and 
has the D^ group symmetry. It consists of hexago- 
nal graphite-like sheets but with an ABAB- • ■ stack- 
ing with boron atoms in layer A found directly be- 
low nitrogen atoms in layer B. The ground state prop- 
erties of h-BN are obtained by means of the all- 
electron Projector- Augmented- Wave method (PAW) [l7j 
from density-functional theory (DFT) in the local-density 
approximation (LDA)^], yielding LDA band-structure 
energies, and wave functions needed to calculate all 
quantities occurring in the GW self-energy operator and 
the electron-hole interaction. The experimental lattice 
parameters 19] (a=4.73 a.u and c/a=2.66) are used to 
perform all calculations. 

The quasiparticle energies are determined by solving 
the QP equation, where the self-energy operator S is 
written down in Hedin's GW approximation |20| which 
includes dynamic polarization in the random-phase ap- 
proximation (RPA). In this case £ can be written as 
the product of the Green's function and the dynami- 
cally screened Coulomb interaction within the RPA. In 
this work, a plasmon pole model|2l| is used to mimic 
the frequency dependence of the dielectric matrix in or- 
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FIG. 1: Calculated electronic band structure along high- 
symmetry lines for bulk hexagonal boron nitride. The thin 
solid lines display the LDA results while the solid lines with 
open circles show the GWA results. The energy scale is rela- 
tive to the top of the valence band maximum (VBM) located 
at the Ti point near K along the F-K high symmetry direc- 
tion. 



der to reduce the computational cost of the GW calcula- 
tions. More details about our implementation are given 
elsewhere [l4|. 

In Fig. both the calculated LDA (thin solid lines) 
and GW (solid lines with open circles) band structures 
are plotted along the high symmetry directions. Within 
the LDA, we found an indirect gap of 4.02 eV between the 
bottom of the conduction band at the M point and the 
top of the valence band near K (T\ point). The minimum 
direct gap is located at the H point and is 4.46 eV. Our 
results are in excellent agreement with other LDA cal- 
culations [Til I22I |23| . The main effect of the self-energy 
correction is to increase the gap from 4.02 eV (LDA) to 
5.95 eV and the minimum direct gap from 4.46 eV (LDA) 
to 6.47 eV. Our GW results are in good agreement with 
the calculations of Capellini et al 13] who found values 
6.04 and 6.66 eV, respectively. Figure Q shows also that 
the GW corrections to the LDA are strongly dependent 
on the degree of orbital localization and that the overall 
corrections can not be reproduced by a rigid band shift 
of the conduction states with respect to the top of the va- 
lence states. Such a behavior has already been reported 
for layered compounds like graphite|24( or LiBC |25|. 

The wide range of experimental valuesQ, 0, H, 0, El 
reported for the band gap (between 3.2 and 5.97 eV) may 
be due, in large part, to the difficulty to synthesize high 
quality h-BN crystals or to the difficulty to extract the 
band gaps from experiment. The most accurate results 
have been obtained by Watanabe et aLpJ who resolved 
the excitonic structure within the gap at very low tem- 
perature (8 K) and deduced from the identification of a 
plateau in their absorption spectrum that the QP gap is 
5.971 eV. Notice that the inferred band gap is very close 
to our calculated indirect band gap. It's worth point- 



ing out that the interpretation of the optical spectra of 
non-conducting solids is very intricate because of strong 
electron- hole correlation effects. 

In the absence of exciton-phonon coupling, emission 
or absorption of a photon involves excited states l^q) of 
total wavevector q corresponding to the photon momen- 
tum, hereafter approximated as q = 0. Due to electron- 
electron interaction, these excited states can be viewed 
as a linear combination of quasielectron-quasihole pairs: 

l*q=0) = E ^ck^klO) (1) 

k,c,t; 

where a^ k creates a quasielectron and a v k a quasihole in 
the many-body ground state |0). The electron-hole am- 
plitudes A^ cW and the excitation energies E x are obtained 
by solving an effective two-particle Schrodinger equat ion, 
which originates from the Bethe-Salpeter eauation|26j 
and provides the key ingredients to calculate the imagi- 
nary part of the dielectric function including both local 
field (LF) and excitonic effects pll27|. 

As the screening controls the strength of the electron- 
hole interaction, we inverted directly the dielectric func- 
tion at the RPA level to obtain the static dielectric func- 
tion £oo(q — > 0). The LDA, without LF effects, produces 
values of 4.67 and 2.98 for E 1 c and E || c, respec- 
tively, where E is the electric field vector. As expected, 
the LF effects push the oscillator strength towards higher 
energies |28| reducing the static dielectric function values 
to 4.40 and 2.53, respectively. Notice that the LF have 
a stronger effect for the out of plane polarization, i.e., 
E || c. This is due to the rather large inhomogeneity 
of the h-BN charge density along the c direction. These 
results show that the screening is rather anisotropic and 
that excitonic effects are expected to be important, i.e., 
the smaller dielectric constant, the larger the excitonic 
effects due to the poor screening of the electron-hole in- 
teraction. 

The calculated optical absorption spectra for h-BN 
with and without excitonic effects are shown in Fig. [21 
and are compared to the experimental results[To| for in- 
plane polarization (E 1 c). The LDA- RPA (dashed line) 
and the GW-RPA (thin solid line) dielectric functions are 
red-shifted and blue-shifted with respect to experiment, 
respectively. Notice that the GW-RPA spectra can not 
be obtained by simply shifting the LDA-RPA spectra to- 
wards higher energies because the corrections to the LDA 
eigenvalues are far from being uniform across the BZ. 
The GW-RPA spectrum for the in-plane polarization ex- 
hibits two prominent structures located at 7.65 eV and 
16.65 eV, respectively. The first structure originates from 
the 7r — * 7r* interband transitions between two parallel 
bands along the M-L high symmetry direction and is de- 
noted by solid arrows in Fig. For such transitions, 
the charge densities for the valence and the conduction 
states both have a p z character and are located on ni- 
trogen and boron atoms, respectively. The second peak 
originates from a — > a* interband transitions. One of 
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FIG. 2: Calculated imaginary part of the dielectric function 
of bulk hexagonal boron nitride with (thick solid line) and 
without (dashed line and thin solid line) excitonic effects for 
light polarized perpendicular (in-plane polarization) and par- 
allel to the c-axis versus photon energy compared with the 
experimental data (empty circles) of Ref . [Io| . The dashed 
lines and the thin solid lines are calculated using the LDA 
and the GW eigenvalues, respectively. 



the corresponding interband transition, occurring at the 
M point, is depicted by a dashed arrow on Fig. ^ The 
GW-RPA spectrum for out-of-plane polarization (E || c) 
displays three major structures located at 12 eV, 12.75 
eV and 16.5 eV and originate from 7r — > a* or a — > 7r* 
interband transitions. The dot-dashed and dotted arrows 
on Fig. n re P resen t one °f the interband transitions con- 
tributing to the peaks located at 12.75 eV and 16.5 eV, 
respectively. 

The assignment of the optical transitions for both po- 
larization directions can be understood by invoking sym- 
metry arguments. Since there is a mirror plane within 
each layer, it's possible to separate bands into states 
of even (a) and odd (it) parity. Therefore, only transi- 
tions between bands of the same symmetry (n — > n* and 
(j — ► (7* ) are allowed for in- plane polarization (E _L c), 
while only transitions between bands of different parity 
(ir — > a* and a — > 7r*) are allowed for out-of-plane polar- 
ization. Notice that this selection rule is only valid for 
k points whose component along the c* axis is strictly 
zero. Thereby, we can explain the rather small but non 
vanishing oscillator strength below 10 eV for the out-of- 
plane polarization. Figure [21 also shows the calculated 
spectra including excitonic effects (thick solid line). The 
electron-hole attraction shifts the oscillator strength to- 
wards lower photon energies by about 1 eV, bringing the 
calculated spectra in better agreement with experiment. 
A bound exciton shows up at 5.85 eV in the calculated 
spectra while a prominent peak occurs at 6.08 eV in the 
experimental spectrum. It is interesting to point out that 
a peak around 6.15 eV, characteristic of h-BN, has been 
found in the absorption spectrum of single-wall boron 




FIG. 3: Two dimensional projections of the probability den- 
sity |* A (r h ,r e )| 2 for the exciton state at 5.780 eV. The left 
panel shows the distribution of the electron relative to the 
hole located slightly above a nitrogen atom (black circle) in 
a plane containing the hole and perpendicular to the c axis. 
The right panel shows the distribution of the electron relative 
to the hole in a plane parallel to the c axis passing through 
the line shown in the left panel. 



nitride nanotubes|29|. A detailed analysis of our calcu- 
lated peak shows that it's made of four excitonic lines 
located at 5.75 , 5.78 , 5.82 and 5.85 eV. Each of this 
excitonic state is built from a coherent superposition of 
quasielectron-quasihole pairs corresponding to n — > it* 
transitions between the two last valence bands and the 
two first conduction bands for k points laying in the 
HLMK plane. The k points which contribute the most to 
these excitonic lines are distributed along directions par- 
allel to the c* axis and in the direction HK, where the 
bands are nearly parallel (see Fig. Notice that these 
excitonic peaks are only visible for in-plane polarization 
because bands of the same parity (ir — > ir*) are involved 
in these excited states. Again, it's worth comparing our 
results with those of Watanabe et a/.Q] who found four 
excitonic lines located at 5.822, 5.945, 5.962 and 5.968 
eV, respectively. These lines were attributed to a se- 
ries of s-like excitons (n = 1, ■ • ■ ,4) by using a Wannier 
model and estimating the band gap to be 5.971 eV. The 
overall agreement with our results is pretty good when 
we shift our calculated excitonic lines by 0.07 eV towards 
higher energy. Nevertheless, we do not reproduce the ex- 
perimental spacing between the different lines. Such a 
difference might be due to the approximations made to 
solve the Bethe-Salpeter equation. Regarding the inter- 
pretation of the excitonic structures inside the gap, our 
results strongly disagree with those of Watanabe et al 
who inferred a binding energy of 0.149 eV from their opti- 
cal absorption spectra. Indeed, our results show that the 
binding energy is much larger, 0.72 eV, and that excitons 
in h-BN belong to the Frenkel and not to the Wannier 
class. 

To address the localization of the excitons for h-BN, 
we calculate the exciton probability |\l/ A (rh, r e )| 2 for a 
given hole position rh and show the distribution of the 
electron relative to the hole for the excited state at 5.78 
eV. Figure |3 displays such a distribution for a hole local- 
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ized slightly above a nitrogen atom. The left panel shows 
that the electron is delocalized on boron atoms and that 
the probability of finding the electron on a boron atom is 
maximum for the nearest neighbors of the nitrogen atom 
where the hole is localized. The right panel shows that 
the electron is confined within the layer where the hole is 
located, which is expected because the screening is rather 
poor along the c direction. The distribution of the elec- 
tron relative to the hole also reflects the orbital character 
of the conduction states which have a p 2 character on B 
atoms. Similar plots have been obtained for other bound 
excitons, confirming that the Frenkel exciton-like picture 
is valid. 

In this letter, results concerning the excited states of 
h-BN have been presented. By calculating the quasipar- 
ticle band structure within the GW approximation, we 
found that h-BN is an indirect gap semiconductor with 
a band gap energy of 5.95 eV in surprisingly close agree- 
ment with the supposed direct gap of 5.971 eV reported 
by Watanabe et aZpJ while our calculated minimum di- 
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